ENTROPIC NONEXTENSIVITY: A POSSIBLE MEASURE OF COMPLEXITY 
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"Beauty is truth, truth beauty,"— that is all 
Ye know on earth, and all ye need to know. 

(John Keats, May 1819) 



An updated review M of nonextensive statistical mechan- 
ics and thermodynamics is colloquially presented. Quite nat- 
urally the possibility emerges for using the value of q — 1 
(entropic nonextensivity) as a simple and efficient manner to 
provide, at least for some classes of systems, some charac- 
terization of the degree of what is currently referred to as 
complexity || . A few historical digressions are included as 
well. 

Keywords: Nonextensive statistical mechanics; En- 
tropy; Complexity; Multifractals; Power laws. 



I. INTRODUCTION 

Would you say that Newtonian mechanics is univer- 
sal? Wc almost hear a big and unanimous yes coming 
from practically all ninctccnth-ccntury physicists. And 
yet, they would all be wrong! Indeed, we know today 
that if the involved masses are very small, say that of 
the electron, Newtonian mechanics has to be replaced by 
quantum mechanics. And if the involved speeds are very 
high, close to that of light in vacuum, special relativ- 
ity has to be used instead. And if the masses are very 
large, like that of the Sun or of a galaxy, general relativ- 
ity enters into the game. It is our present understanding 
that only in the (h, 1/c, G) -> (0, 0, 0) limit, Newto- 
nian mechanics appears to be strictly correct. Now, what 
about Boltzmann-Gibbs statistical mechanics and stan- 
dard thermodynamics? Are they universal? An attempt 
to answer such question and to clarify its implications 
constitutes the main goal of the present effort. Let us 
see what Ludwig Boltzmann has to tell us about this. 
In the first page of the second part of his Vorlesungen 
uber Gastheorie, he qualifies the concept of ideal gas by 
writing: " When the distance at which two gas molecules 
interact with each other noticeably is vanishingly small 
relative to the average distance between a molecule and 
its nearest neighbor - or, as one can also say, when the 
space occupied by the molecules ( or their spheres of ac- 
tion) is negligible compared to the space filled by the gas - 
Were Boltzmann our contemporary, he would have 
perhaps told us that he was addressing systems involving 
what we nowadays call short range interactions! Several 
decades ago, Laszlo Tisza H, in his Generalized Ther- 



modynamics, writes "The situation is different for the 
additivity postulate [...], the validity of which cannot be 
inferred from general principles. We have to require that 
the interaction energy between thermodynamic systems be 
negligible. This assumption is closely related to the homo- 
geneity postulate [. . .}. From the molecular point of view, 
additivity and homogeneity can be expected to be reason- 
able approximations for systems containing many parti- 
cles, provided that the intramolecular forces have a short 
range character" , when referring to the usual thermo- 
dynamic functions such as internal energy, entropy, free 
energy, and others. Also, Peter T. Landsberg [Q, in his 
Statistical Mechanics and Thermodynamics, writes "The 
presence of long-range forces causes important amend- 
ments to thermodynamics, some of which are not fully 
investigated as yet. ". More, during a visit to Rio de 
Janeiro some years ago, he told me that such restric- 
tions should be in the first pages of all books on thermo- 
dynamics... but they are not! The title of the 1993 Na- 
ture article Q by John Maddox focusing on black holes 
is "When entropy does not seem extensive". What is 
the viewpoint on such matters of our colleagues math- 
ematicians? Distinguished expert in nonlinear dynami- 
cal systems, Floris Takens writes || "The values of pi 
are determined by the following dogma: if the energy 
of the system in the i th state is Ei and if the temper- 
ature of the system is T then: pi = exjp{—Ei/kT}/Z(T), 
where Z(T) = YJ f exp{— Ei/kT} [...] This choice of Pi is 
called Gibbs distribution. We shall give no justification 
for this dogma; even a physicist like Ruelle disposes of 
this question as " deep and incompletely clarified". " It 
is known that mathematicians sometimes use the word 
dogma when they do not have the theorem! As early as 
1964, Michael E. Fisher , and later on with David Ru- 
elle as well as with Joel L. Lebowitz, addressed in detail 
such questions They established, for instance, that 
a d-dimensional classical system (say a fluid) including 
two-body interactions that are nonsingular at the origin 
and decay, at long distances, like r~ a , exhibits standard 
thermostatistical behavior if a/d > 1. An interesting 
question remains open: what happens when a/d < 1 
? This is the case of all self-gravitating systems (e.g., 
galaxies, black holes), a fact which explains well known 
difficulties in using standard thermodynamics to address 
them. These difficulties are not without relation with 
the case of a single hydrogen atom (d = 3 and a = 1), 
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for which no tractable Boltzmann-Gibbs thermostatis- 
tical calculations are possible (diverging partition func- 
tion) unless we confine the hydrogen atom in a box ||. 
Of course, no such difficulties are encountered if we fo- 
cus on a single harmonic oscillator, or on a single spin in 
an external magnetic field, or on a Lennard- Jones fluid 
(d = 3 and a = 6), or even on a neutral plasma, where 
the Coulombian interaction is "dressed" thus becoming 
a thermostatistically innocuous, exponentially decaying 
interaction. No severe anomalies appear in such cases, 
all relevant sums and integrals being finite and perfectly 
defined. More generally speaking, a variety of recent 
analytical and numerical results (see || and references 
therein) in classical systems (with no singularities at the 
potential origin or elsewhere) suggest that, at the appro- 
priately defined ultimate t — > oo thermal equilibrium, the 
internal energy per particle grows, when the number of 
particles N grows, like TV = [JV 1 - a / d - a/d]/[\ - a/d]. 
We easily verify that N approaches a finite constant if 
a/d > 1, diverges logarithmically if a/d — 1, and di- 
verges like iV 1_Q / d if < a/d < 1, thus recovering, for 
a = 0, the well known Molecular Field Approximation 
scaling with N (usually englobed within the coupling 
constant whenever dealing with magnetic models such 
as the Ising ones, among many others). We trivially see 
that the standard thermodynamical extensivity is lost in 
such cases. The physical behavior for the marginal case 
a/d = 1 is intimately related to the mathematical fact 
that some relevant sums that are absolutely convergent on 
the extensive side and divergent on the nonextensive one, 
become conditionally convergent in that case. This en- 
ables us to understand why the amount of calories to be 
provided to a table in order to increase its temperature 
in one degree only depends on its weight and material 
(iron, wood), whereas the amount of Coulombs we must 
provide to a capacitor to generate a one Volt potential 
difference at its ends also depends on its shapel Indeed, 
the relevant interaction in the latter being the perma- 
nent dipole-dipole one (hence d = a = 3), the capacity 
depends on whether the capacitor is a slab, round, cylin- 
drical or otherwise. Such behavior dramatically contrasts 
with the extensivity usually focused on in textbooks of 
thermodynamics. In fact, it is after all quite natural that 
self-gravitating systems provided the first physical appli- 
cation of the formalism we are addressing here. Indeed, 
physicist Angel Plastino and his son astronomer Angel R. 
Plastino showed in 1993 |t0| that by sufficiently depart- 
ing from BG thermostatistics, it is possible to overcome 
an old difficulty, namely to have the physically desirable 
feature that total energy, entropy and mass be simulta- 
neously finite. 

Let us now try to deepen our quest for understanding 
the restrictions within which Boltzmann-Gibbs statistical 
mechanics is the appropriate formalism for describing a 
(fantastic!) variety of systems in nature. Are long-range 
interactions the only mechanism which creates anoma- 
lies? The answer seems to be no. Indeed, several indica- 



tions exist which suggest that mesoscopic dissipative sys- 
tems (like granular matter) also are thermostatistically 
anomalous: the velocity distribution measured in a va- 
riety of computational experiments is, as argued by Leo 
P. Kadanoff et al 0, Y.-H. Taguchi and H. Takayasu 
IP^I , Hans J. Herrmann 13 and others, far from be- 
ing Maxwell's Gaussian. On general grounds, multifrac- 
tally structured systems also tend to exhibit anomalies. 
And I believe that it would be surprise for very few that 
the same happened with systems based on strongly non- 
markovian microscopic memory (e.g., a memory function 
decreasing towards the past as a slow enough power-law) . 

As an attempt to theoretically handle at least some of 
these anomalies, I proposed in 1988 (lj] the generaliza- 
tion of the Boltzmann-Gibbs formalism by postulating a 
nonextensive entropy S q (See Sections I and II) which 
recovers the usual logarithmic measure as the q = 1 par- 
ticular case. Before going further, let us describe S q , its 
properties and historical context. 



II. THE NONEXTENSIVE ENTROPY Sq 

A. Properties 

In the case of a discrete number W of microscopic con- 
figurations, S q is given by [Q 
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(similar expressions correspond to the cases where we 
have a continuum of microscopic configurations, or when 
the system is a quantum one). 



In the limit q — > 1 we have p. 



In pi 



l)lnpi, hence Si — —k^2i—iPi In Pi, which is the usual 
Boltzmann-Gibbs-Shannon expression. Also, it can be 
shown that 

(i) S q > 0, and, for q > 0, equals zero in the case 
of certainty (i.e., when all probabilities vanish but one 
which equals unity); 

(ii) S q attains its extremum (maximum for q > and 
minimum for q < 0) for equiprobability (i.e., pi — l/W), 
thus becoming S q — klW 1 ^' 1 — l]/[l — q\. This expression 
becomes, for q = 1, Si = khiW, which is the formula 
graved on Boltzmann's tombstone in the Central Ceme- 
tery in Vienna; 

(iii) S q is concave (convex) for q > (q < 0), which is 
the basis for thermodynamic stability; 

(iv) If two systems A and B are independent in the 
sense of the theory of probabilities (i.e., pf^ 3 — pf x pf), 
then 

S q (A + B) 

k 



S q (A) 



S g (B) 



+ (1-?) 



Sq(A) Sq(B) 



(2) 
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hence, superextensivity, extensivity or subextensivity oc- 
curs when q<l, q = lorq>l respectively; 

(v) If the W states of a system are divided into Wl 
and Wm states (Wl + Wm = W), and we define pl — 

J2 WL statesPl and VM = E Wm states P™ (PL + PM = 1), 

then 

Sq({pi}) = S {pl,Pm) + p q L S q {{pi/p L }) + p q M S q ({p m /p M 

(3) 

which generalizes Shannon's celebrated formula; 

(vi) It can be easily checked that — d(J2iPi )/dx\ x =i = 
— pj In pj. Sumiyoshi Abe, in Funabashi, made in 
1997 |15|] a simple but very deep remark. He considered 
the so-called Jackson derivative. The mathematician F. 
Jackson generalized, one century ago, the usual differen- 
tial operator. He defined D q f(x) = [f(qx) — f(x)]/[qx — 
q], which reproduces d/dx in the limit q — ► 1. Abe' s ob- 
servation consists in the following easily verifiable prop- 
erty, namely — D q (Y^ i pf )/dx\ x =i — S q ! Since it is clear 
that Jackson derivative "tests" how the function f(x) 
"reacts" under dilatation of x in very much the same 
way the standard derivative "tests" it under translation 
of x, Abe's remark opens a wide door to physical insight 
onto thermodynamics. It was so perceived, I believe, by 
Murray Gell-Mann, of the Santa Fe Institute. Indeed, 
one year ago, during a conference in Italy, he asked me 
whether I would go to conjecture that systems with rele- 
vant symmetries other than say translation or dilatation 
invariances would need entropies other than Si or S q . 
My answer was probably yes. Indeed, there are symme- 
tries that appear in general relativity or in string theory 
that are characterized by a sensible amount of parame- 
ters. The concept of thermodynamic information upon 
such systems, i.e., the appropriate entropy could well be 
a form involving essentially a similar number of param- 
eters. If so, then symmetry would control entropy, and 
S q would only be the beginning of a presumably long 
and fascinating story! Since we have learnt from Gell- 
Mann and others, how deeply symmetry controls energy, 
we could then say that symmetry controls thermodynam- 
ics, that science which nothing is but a delicate balance 
between energy end entropy. More precisely, symmetry 
would then determine (see also p6| ) the specific micro- 
scopic form that Clausius thermodynamic entropy would 
take for specific classes of systems. Symmetry controlling 
thermodynamics! If such a conjectural scenario turned 
out to be correct, this would not be to displease Plato 
with his unified view of truth and beauty! 

(vii) Along the lines of Claude E. Shannon's theorem, 
it has been shown by Roberto J. V. dos Santos, in Ma- 
ceio, Brazil [JL7J, that a set of conditions generalizing 
(through properties (iv) and (v)) that imposed by Shan- 
non uniquely determines S q . Analogously, Abe recently 
showed [|18[ that consistently generalizing Khinchin's set 
of conditions, once again a unique solution emerges, 
namely S q . 



B. Brief review of the labyrinthic history of the 
entropies 

The entropy and the entropic forms have haunted 
physicists, chemists, mathematicians, engineers and oth- 
ers since long! At least since Clausius coined the word 
and wrote "The energy of the universe is constant, its 
entropy tends to increase" (See Appendix). The whole 
story is a labyrinth plenty of rediscoveries. After Ludwig 
Boltzmann and Josiah Willard Gibbs introduced and first 
analyzed, more than one century ago, their respective ex- 
pressions |^(| of the entropy in terms of microscopic quan- 
tities, many generalizations and reformulations have been 
proposed and used. Let us mention some of the most rel- 
evant ones. John von Neumann introduced a quantum 
expression for the entropy in terms of the density opera- 
tor; when the operator is diagonalized, the traditional ex- 
pression (herein referred, because of the functional form, 
to as the Boltzmann-Gibbs entropy, although these two 
scientists used to work with different, though consistent, 
expressions is recovered. Shannon reinvented, more 
than half a century ago, a binary version of the B-G func- 
tional form and interpreted it in terms of information and 
communication theory. It is no doubt very interesting to 
notice that in The mathematical theory of communication 
he wrote, in reference to the specific form he was using, 
that " This theorem, and the assumptions required for its 
proof, are in no way necessary for the present theory. It 
is given chiefly to lend a certain plausibility to some of 
our later definitions. The real justification of these def- 
initions, however, will reside in their implications. ". In 
1957, Edwin Jaynes attributed to the whole formalism a 
very generic, and still controversial in spite of an enor- 
mous amount of practical applications, anthropomorphic 
shape. Since those years, people working in cybernetics, 
information theory, complexity and nonlinear dynamical 
systems, among others, have introduced close to 20 dif- 
ferent entropic forms! This number has almost doubled, 
this time mainly because of physicists, since the 1988 
paper [Q appeared. 

In 1960, A. Renyi proposed [^lj a form which recovers 
that of Shannon as a special case. His form is always ex- 
tensive, but not always concave (or convex) with regard 
to the probability distributions. It seems that it was a 
rediscovery; indeed, according to I. Csiszar's 1974 critical 
survey p2| , that form had essentially already been intro- 
duced by Paul-Marcel Schutzenberger in 1954 |^3). In 
1967, presumably for cybernetic goals, J. Harvda and F. 
Charvat introduced p4| , although with a different prefac- 
tor (fitted for binary variables) , the entropic form herein 
noted S q , and studied some of its mathematical prop- 
erties. One year later, I. Vajda further studied plflthis 
form, quoting Harvda and Charvat. In 1970, in the con- 
text of information and control, Z. Daroczy rediscovered 
po| this form (he quotes neither Harvda and Charvat, nor 
Vajda). Later on, in 1975, B.D. Sharma and D.P. Mittal 
introduced pTJ and studied some mathematical proper- 
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ties of a two-parameter form which recovers both Renyi's 
and S q as particular instances! In his 1978 paper in Re- 
views of Modern Physics, A. Wehrl mentions [p8[ the 
form of S q , quoting Daroczy, but not Harvda, Charvat, 
Vajda, Sharma and Mittal. In 1988, in the Physics com- 
munity, a new rediscovery, by myself this time. Indeed, 
inspired by mult ifr act als, I postulated |14| the form S q 
as a possible new path for generalizing Boltzmann-Gibbs 
statistical mechanics; through optimization of S q I ob- 
tained an equilibrium-like distribution which generically 
is a power-law, and reproduces Boltzmann's celebrated 
factor as the particular limit q — > 1. In my 1988 paper, 
I quoted only Renyi's entropy, having never taken notice 
of all these generalized forms, whose existence (and that 
of many others!) I have been discovering along the years. 
In fact, Renyi's entropy can be shown |3 to be a mono- 
tonically increasing function of S q . However, concavity is 
not preserved through monotonicity. So what?, since the 
optimizing distribution will be exactly the same. Well, it 
happens that statistical mechanics is much more than an 
equilibrium-like distribution optimizing an entropy under 
some generic constraints! It includes a variety of other 
relevant points, such as the role played by fluctuations, 
linear responses, thermodynamical stability, to name but 
a few. Nobody would, I believe, try to reformulate BG 
statistical mechanics using, instead of the BG entropy, 
say its cube, just because it is a monotonically increasing 
function of the BG one! Let us also mention that, within 
various scientific and technological communities, close to 
30 (!) different entropic forms have been considered along 
the last half a century. 

Let us summarize some of what we have said. The 
BG entropy, Si, is both concave and extensive, the Renyi 
entropy, Sff, is extensive but nonconcave (for all q > 0; 
concave only for < q < 1), the one focused on herein, 
S q , is nonextensive but concave (for all q > 0), and the 
Sharma-Mittal entropy is nonextensive and nonconcave. 
The latter is determined by two parameters and contains 
the other three as particular cases. The Renyi entropy 
and S q are determined by one parameter (q) and share a 
common case, namely the BG entropy, which is nonpara- 
metric. Since Renyi's, Sharma-Mittal's and the present 
entropy depend from the probability set through one and 
the same expression, namely they are all three re- 

lated through simple functions. In particular, 



q l-o 



= k 



ln[l + (1 - q)S q /k} 



1 



(V?), (4) 



with = Si = -kJ2iPi^ n Pi- 

For completeness, let us close this section by men- 
tioning that expressions of the type J^iPi have also 
been, long ago, discussed on mathematical grounds by 
Hardy, Littlewood and Polya [|l9| (and, for q — 2, by the 
Pythagoreans!). 



III. EQUILIBRIUM-LIKE STATISTICS 

The next step is to follow along Gibbs path. For in- 
stance, to formally obtain the "equilibrium" distribution 
associated with the canonical ensemble (i.e., a conserva- 
tive system at "equilibrium" with a thermostat, where 
the physical meaning of this state will turn out to be, as 
we shall see later on, close to a stationary one), we need 
to impose a constraint on total energy. We shall adopt 
p9| that (H) q = J2i^i e i = where U q is fixed and 
finite, {cj} is the set of energy levels of the system Hamil- 
tonian Ji, and the escort distribution {Pi} is defined as 
Pi = pf /Ej Pj]) from now on we refer to {...) q as the 
normalized q-expectation value. At the present stage, a 
plethora of mathematical reasons exist for adopting this 
particular generalization of the familiar concept of inter- 
nal energy. Not the least of these reasons is the fact that 
such definition enables relevant sums and integrals to be 
finite, which would otherwise diverge (e.g., the second 
moment of Levy distributions diverges, whereas the sec- 
ond q- moment with q appropriately chosen is finite). We 
can say that in this sense the theory becomes regularized. 
The optimization of S q with this constraint (to which we 
associate a Lagrange parameter (3) leads to 

Pi = [1 - (1 - q)^ - Uqf^/Zq (/}, EE 

(5) 

with Z q ee £.[1 - (1 - q)(3q( ej - tf,)]VCi-d (See Fig. 
1). We easily verify that q — ► 1 recovers the celebrated, 
exponential Boltzmann factor 



/Zi (Z 1 = 



)• 



(6) 



For q > 1 a power-law tail emerges; for q < 1 the formal- 
ism imposes a high-energy cutoff, i.e., pi = whenever 
the argument of the power function becomes negative. 
One comment is worthy: this distribution is generically 
not an exponential law, i.e, it is generically not factor- 
izable (under sum in the argument), and nevertheless is 
invariant under the choice of zero energy for the energy 
spectrum! (this is in fact one of the aside benefits of 
defining the constraints in terms of normalized distribu- 
tions like the escort ones). 

Eq. (5) can be rewritten in the following convenient 
form: 



Pi 



[l-(l-q)P'e l ] 1 /^/Z q 



(7) 



with Z' q 



11/(1-9) 



where /3' is a simple 



Ej[i-(i-«)^j 

function of [3 q and Up- 

It can be shown [B9[, along the lines of the 1991 pa- 
per with Evaldo M.F. Curado ^9|, that this formalism 
satisfies, for arbitrary q, 1/T = dS q /dU q (T = l/(k/3)), 
and the entire Legendre transform structure of standard 
thermodynamics. Moreover, the relations which provide 
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U q and F q = U q — TS q in terms of Z q (i.e., the con- 
nection between the microscopic and macroscopic de- 
scriptions of the system), remain basically the same as 
the usual ones, the role played by the logarithmic func- 
tion In a: being now played by [x 1-9 — 1)/(1 — q)- Fi- 
nally, various important thcrmostatistical theorems can 
be shown to be q-invariant. Among them, we must men- 
tion the Boltzmann H-theorem (macroscopic time irre- 
versibility, i.e., essentially the second principle of ther- 
modynamics, first tackled within this context by Ana- 
nias M. Mariz in Natal, Brazil J|(|), Ehrenfest theorem 
(correspondence principle between classical and quantum 
mechanics), Einstein 1910 factorizability of the likelihood 
function in terms of the entropy (this q-generalization 
was first realized in 1993 during a private discussion 
with Manuel O. Caceres, from Bariloche, Argentina), On- 
sager reciprocity theorem (microscopic time reversibil- 
ity), Kramers and Wannier relations (causality), and, 
apparently, even Pesin relation for nonlinear dynamical 
systems, although the last one remains at the level of 
plausibility (no proof available). The q- generalization of 
the fluctuation-dissipation theorem was first addressed 
by Anna Chame and Evandro V.L. de Mello, in Niteroi, 
Brazil pi] ], and then, along with the generalization of 
Kubo's linear response theory, by A.K. Rajagopal, in the 
Naval Research Laboratory in Washington [32]. Since 
then, a variety of important statistical mechanical ap- 
proaches such as the Bogolyubov inequality, the many- 
body Green function and path integral formalisms have 
been generalized by Renio S. Mendes in Maringa, Brazil, 
Ervin K. Lenzi in Rio de Janeiro, and collaborators fl32| . 
Nonextensive versions of Fermi-Dirac and Bose-Einstein 
statistics have been discussed [|3| by Fevzi Buyukkilic, 
Dogan Demirhan and Ugur Tirnakli in Izmir, Turkey, 
and their collaborators. Finally, it is certainly impor- 
tant to mention that, very recently, Rajagopal and Abe 
have done an exhaustive review [|34f| of apparently all 
the traditional methods for obtaining the BG equilibrium 
distribution, namely the Balian-Balazs counting method 
within a microcanonical basis, the Darwin-Fowler steep- 
est descent method, and the Khinchin's large-numbers 
method. They g-generalized all three, and systemati- 
cally obtained the same power-law distribution obtained 
fl4] , p9"[ above through the Gibbs' method, i.e., the opti- 
mization of an appropriately constrained entropy. 

In order to see the nonextensive formalism at work, 
let us address Q an important question which long re- 
mained without satisfactory solution, namely, why Levy 
distributions are ubiquitous in nature, in a manner which 
by all means is similar to the ubiquity of Gaussians? 
To introduce the point, let us first address the Gaussian 
case. If we optimize the BG entropy Si by imposing a 
fixed and finite average (x 2 ), we obtain a Gaussian one- 
jump distribution. By convoluting N such distributions 
we obtain a A-jump distribution which also is a Gaus- 
sian. From this immediately follows that (x 2 )(N) oc N. 
If we consider that N is proportional to time, we have 
the central result of Einstein's 1905 celebrated paper on 



Brownian motion, which provided at the time strong ev- 
idence in favor of Boltzmann's ideas. We may reword 
this result by saying that the foundation of Gaussians in 
nature lies upon two pillars, namely the BG entropy and 
the central limit theorem. What about Levy distribu- 
tions? The difficulty comes from the fact that all Levy 
distributions have infinite second moment, due to their 
long tail (oc l/|a;| 1+7 with < 7 < 2). Consequently, 
what simple auxiliary constraint to use along with the 
optimization of the entropy? In 1994, Pablo Alemany 
and Damian H. Zanette in Bariloche, Argentina, settled 
p6| the basis for the solution, namely optimization of S q 
while imposing a finite g-expectation value of x 2 . The 
one-jump distribution thus obtained is proportional to 
[1 — (1 — q)f3x 2 } 1 /( 1 ~ q \ which respectively reproduces the 
Gaussian, Lorentzian (Cauchy for the mathematicians) 
and completely flat distributions for q = 1, 2 and 3. For 
q < 1 the support is compact, and for q > 1 a power- 
law tail is obtained. The proof was completed one year 
later |$5|, when it was argued that the Levy-Gnedenko 
central limit theorem was applicable for q > 5/3, thus 
approaching, for N » 1 precisely Levy distributions for 
the A-jump distribution. We may then summarize this 
result by saying that the foundation of Levy distributions 
in nature also has two pillars, namely S q and the Levy- 
Gnedenko theorem. It follows that 7 = 2 for q < 5/3 and 
7 = (3 — q)/{q — 1) for 5/3 < q < 3. Since 7 can also 
be interpreted as the fractal dimension to be associated 
with Levy walks, this turned out to be the first exact con- 
nection between nonextensive statistical mechanics and 
scaling. Similar ideas have been successfully put forward 
by Plastino and Plastino |37j and by Lisa Borland |3|] 
and others, concerning nonlinear as well as fractional- 
derivative Fokker-Planck-like equations. In all these var- 
ious types of anomalous diffusion, the Gaussian solution 
associated with Jean Baptiste Joseph Fourier's bicenten- 
nial heat equation is recovered as the q = 1 limit. Let 
me also mention that Hermann Haken, in Stuttgart, and 
collaborators have recently applied related ideas within 
a theoretical model of the human brain |3S| ]. 

In the same vein which led to the above discussion of 
Levy distributions, if we optimize S q imposing, besides 
normalization of p{x), finite values for the normalized 
q-expectation values of both x and x 2 , we obtain 

p q (x) oc [1 - (1 - q)((3 {1) x + p^x 2 )] 1 '^ , (8) 

where (3^ and /3' 2 ' are determined by the constraints. 
This distribution contains, as particular cases, a consid- 
erable amount of well (and not so well) known distribu- 
tions, such as the exponential, the Gaussian, the Cauchy- 
Lorentz, the Edgeworth [fo|](for /jM = 1, (3^ = 0, and 
x > 0), the r- and the Student's s- ones jyj, among 
others. Let us end with a short and amusing historical 
remark concerning that which can be considered as the 
most famous of all nontrivial distributions, namely the 
normal one. In contrast with what almost everybody 
would naively think, it was |l2| first introduced by Abra- 
ham De Moivre in 1733, then by Pierre Simon de Laplace 
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in 1774, then by Robert Adrain in 1808, and finally by 
Carl Friedrich Gauss in 1809, nothing less than 76 years 
after its first publication! 

IV. MICROSCOPIC DYNAMICAL FOUNDATION 

To have a deeper understanding of what this general- 
ized formalism means, it is highly convenient to focus on 
the mixing properties. 

A. Mixing in one-dimensional dissipative maps 

Let us now address a very simple type of system, which 
will nevertheless provide important hints about the phys- 
ical and mathematical significance of the present formal- 
ism. Let us consider the well known logistic map, namely 

x t+ i = l~ax 2 t (t = 0, 1, 2, < a < 2) (9) 

For a < a c = 1.40115519..., the attractors are finite cy- 
cles (fixed points, cycle-2, cycle-4, etc) and the Lyapunov 
exponent (hereafter noted Ai) is almost everywhere neg- 
ative (i.e., strongly insensitive to the initial conditions 
and any intermediate rounding). For a c < a < 2, the at- 
tractor is chaotic for most values (i.e., it has an infinite 
number of elements, analogously to an irrational num- 
ber) and the Ai exponent is consistently positive (i.e., 
strongly sensitive to the initial conditions and any in- 
termediate rounding). There is however, in the inter- 
val [0, 2], an infinite number of values of a for which 
Ai = 0. What happens with the sensitivity to the ini- 
tial conditions in those marginal cases? Let us focus on 
this problem by first recalling the definition of Ai. If two 
initial conditions Xq and x' Q differ by the small amount 
Aa;o , we can follow the time evolution of Axt through the 
quantity £ = \imAx ^o{Axt/ Axq). The most frequent 
situation is that £ = e Al which defines the Lyapunov 
exponent Ai. We can trivially check that £(t) is the so- 
lution of d^/dt — Ai £. We have conjectured that, 
whenever Ai vanishes for this map (and similar ones), 
the controlling equation becomes d^/dt = X q £ 9 , hence 
£ = [1 + (1 — q)X q t] 1 ^ 1 ^^. This result recovers the usual 
one at the q = 1 limit; also, it defines a generalized Lya- 
punov coefficient X q , which, like Ai, inversely scales with 
time, but now within a law which generically is a power, 
instead of the standard exponential. The two paradig- 
matic cases are q < 1 (with X q > 0), hereafter referred 
to as weakly sensitive to the initial conditions, and q > 1 
(with X q < 0), hereafter referred to as weakly insensitive 
to the initial conditions. Within this unified scenario, 
we would have strong chaos for q — 1 and Ai > (i.e., 
exponential mixing) and weak chaos for q < 1, Ai = 
and A g > (i.e., power-law mixing). All these various 
possibilities have indeed been observed for the logistic 
map. If we start with xq ~ 0, we have (i) Ai < 
for a = a c - 10~ 3 , (ii) A x > for a = a c + 10~ 3 , (iii) 



Ai = 0, q > 1 and X q < for a doubling-period bifurca- 
tion (e.g., a — 3/4) as well as for a tangent bifurcation 
(e.g., a — 7/4), and finally (iv) Ai = 0, q < 1 and X q > 
at the edge of chaos (i.e., a — a c ). The last case is, 
for the present purposes, the most interesting one. This 
power-law mixing was first observed and analyzed by P. 
Grassberger, A. Politi, H. Mori and collaborators p[ . 
Its importance seems to come mainly from the fact that, 
in the view of many authors, "life appeared at the edge 
of chaos"! Moreover, this type of situation seems to be 
a paradigm for vast classes of the so-called complex sys- 
tems, much in vogue nowadays. Let us say more about 
this interesting point. The precise power-law function 
£ indicated above appears to be the upper bound of a 
complex time-dependence of the sensitivity to the initial 
conditions, which includes considerable and ever lasting 
fluctuations. From a log-log representation we can nu- 
merically get the slope and, since this is to be identified 
with 1/(1 — q), we obtain the corresponding value q*. 
For the logistic map we have obtained Q q* — 0.24.... 
If we apply these concepts to a more general map, like 
x t+ i = 1 — a\x t \ z (0 < a < 2; z > 1), we can check, 
and this is long known, that the edge of chaos depends 
on z. More precisely, when z increases from 1 to infinity, 
a c increases from 1 to 2. If for a c we check the sensitiv- 
ity to the initial conditions, we obtain that q* increases 
from — oo to close to 1 (though smaller than 1), when z 
increases from 1 to infinity. Simultaneously, the fractal 
dimension df of the chaotic attractor increases from 
to close to 1 (though smaller than 1). If instead of the 
logistic family of maps, we use some specific family of 
circle maps also characterized by an inflexion exponent 
z, it can be shown that, once again, q* increases with z 
for all studied values of z (z > 3). The interesting fea- 
ture in this case is that df = 1 for all these values of z. 
What we learn from this is that the index q characterizes 
essentially the speed of the mixing, and only indirectly 
how "filled" is the phase space within which this mixing 
is occuring. 

Let us now address, for maps like the logistic one, a 
second manner of obtaining q* , this time based on the 
multifractal geometry of the chaotic attractor. This ge- 
ometry can be, and usually is, characterized by the so 
called f(a) function, a down- ward parabola- like curve 
which generically is below (and tangential to) the bissec- 
tor in the (a, /) space, is concave, its maximal value is df, 
and, in most of the cases, vanishes at two points, namely 
a mm and a max (typically, < a min < d f < a max ). For 
a fractal, one expects, a min = df = a max = f(a rnin ) = 
f(df) — f(ctmax)- For a so-called fat (multi)fractal, df 
attains the value of the Euclidean space within which the 
(mufti) fractal is embedded (this is the case of the fam- 
ily of circle maps mentioned above). Marcelo L. Lyra, 
in Maceio, and myself used (imitating the successful Ben 
Widom style!) in 1998 some simple scaling argu- 
ments and obtained the following relation: 
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This is a kind of fascinating relation since it connects 
the power-law sensitivity of nonlinear dynamical systems 
with purely geometrical quantities. It has been numeri- 
cally checked for at least three different one-dimensional 
dissipative maps. Its analysis for two- or more- dimen- 
sional systems, as well as for conservative ones, would 
of course be very enlightening and remains to be done 
(as do quantum analogs of this problem). Furthermore, 
it does not contain in a generic way the q = 1 limit. 
Indeed, it seems that, for its validity, it is implicitly 
assumed that f(a min ) = f{a max ) — 0. A more gen- 
eral expression could be something like 1/(1 — q*) = 
l/[a min - f(a min )] - l/[a max - f(a max )}; the generic 
q = 1 limit could then correspond to 
The second manner we referred above goes as follows. 
We can numerically construct the f(a) function asso- 
ciated with the chaotic attractor, measure a m i n and 
c*max, and, through the above connection, calculate q*. 
This procedure yields, for the standard logistic map, 
q* = 0.2445..., thus reproducing the number we previ- 
ously obtained from the sensitivity to the initial condi- 
tions. The higher precision we indicate here comes from 
the fact that, for the logistic-like family maps, it is long 
known that a max — z a m i n — > where ap denotes 
the corresponding Feigenbaum constant. This constant 
being known with many digits for z = 2, we can infer q* 
with quite high precision. 

Let us finally address now a third manner of obtaining 
q* , this time directly from the nonextensive entropy S q . 
It is known that, for large classes of nonlinear dynam- 
ical systems, the rate of loss of information as time in- 
creases can be characterized by the so called Kolmogorov- 
Sinai entropy (hereafter noted K\). This quantity is de- 
fined, strictly speaking, through a single-trajectory proce- 
dure, based on the frequencies of appearance, in increas- 
ingly long strips, of symbolic sequences associated with 
partitions of the phase space. However, in apparently 
very large (if not all) classes of problems, this single- 
trajectory procedure can be replaced by an ensemble- 
based one, whose numerical implementation is by far sim- 
pler computationally than the implementation of algo- 
rithms based on sequences. We shall restrict ourselves to 
this ensemble-based procedure. Let us illustrate, on the 
logistic map, the steps to be followed in order to calcu- 
late K\. We choose a value of a for which Ai is positive, 
say a — 2 (whose corresponding Lyapunov exponent is 
Ai = In 2). We then partition the accessible phase space 
(— 1 < x < 1, in our example) inW » 1 little windows, 
and choose (in any way), inside one of those windows, 
TV initial values for xq. We then apply the map onto 
each of those initial conditions, which, as time evolves, 
will start spreading within the accessible phase space. At 
time t, we will observe a set {Ni(t)} of points inside the 

W windows {Y^Li N i( l ) = N > V *)- We next define a set 
of probabilities pi = Ni/N (i — 1, 2, W), which en- 



able in turn the calculation of Si(t) = — J2 i=1 Pi Inpi. 
This entropy starts, by construction, from zero at t = 0, 
and tends to linearly increase as time goes on, finally 
saturating at a value which, of course, cannot exceed 
\nW. The Kolmogorov-Sinai-like entropy rate is then 
defined through K\ = lim^oo limvi/^oo lirnAr— >oo S\(t)/t. 
We expect, according to the Pesin equality, K\ = Ai. In- 
deed, we numerically obtain K\ = In 2 (see Fig. 2(a)). 
We can analogously define a generalized quantity, namely 
K q = Iimt_ >00 limvi/^oc lim— ><» S q (t)/t. And we can nu- 
merically verify that, for say a = 2, K q vanishes (di- 
verges) for any value of q > 1 (q < 1), being finite only 
for q = 1. Let us now turn onto the edge of chaos, i.e., 
a = a c . Vito Latora and Andrea Rapisarda, in Catania, 
Italy, Michel Baranger, at MIT, and myself |4Sfl , have nu- 
merically verified the same structure, excepting for the 
fact that now K q vanishes (diverges) for any value q > q* 
(q < q*), being finite only for q = q* = 0.24... ! (see Fig. 
2(b)). Tirnakli et al 0] have also verified the validity of 
this structure for various typical values of z within the 
logistic- like family of maps. 

Let us summarize this section by emphasizing that we 
have seen that it is possible to arrive to one and the same 
nontrivial value of q* through three different and sugges- 
tive roads, namely sensitivity to the initial conditions, 
multifractal structure, and rate of loss of information. By 
so doing, we obtain insight onto the mathematical and 
physical meaning of q and its associated nonextensive 
formalism. Moreover, a pleasant unification is obtained 
with the concepts already known for standard, exponen- 
tial mixing. A very interesting question remains open 
though, namely the generalization of Pesin equality. It 
might well be that, as conjectured in p3[ , K q equals A g 
(or, more generically, an appropriate average of X q (x), 
whenever nonuniformity is present). However, a math- 
ematical basis for properties such as this one is heavily 
needed. Some recent numerical results along these lines 
by Paolo Grigolini, in Denton, Texas, and collaborators 
P6| seem promising. In the same vein, it is known that 
whenever escape exists, as time evolves, out of a specific 
region of the phase space of a chaotic system, K\ be- 
comes larger than Ai, the difference being precisely the 
escape rate p7| . All these quantities are defined through 
exponential time-dependences. What happens if all three 
are power laws instead? A natural speculation emerges, 
namely that perhaps a similar relation exists for q ^ 1. 

B. Mixing in many-body dissipative maps 

In the previous section we described a variety of 
interesting nonextensive connections occurring in one- 
dimensional dissipative systems. A natural subsequent 
question then is what happens when the dissipative sys- 
tem has many, say N » 1, degrees of freedom? Such 
is the case of essentially all the models exhibiting Per 
Bak's self-organized criticality (SOC). Francisco Tamarit 
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and Sergio Carinas, in Cordoba, Argentina, and myself 
have discussed Qjthis issue for the Bak-Sneppen model 
Q for biological evolution. By numerically studying the 
time evolution of the Hamming distance for increasingly 
large systems, a power-law was exhibited. By identify- 
ing the exponent with 1/(1 — q*), the value q* ~ —2.1 
was obtained. Totally analogous results are obtained for 
the Suzuki-Kaneko model for imitation games (e.g., bird 
singings) |3(|. It would be of appreciable importance if 
the same numbers were achieved through the other two 
procedures, namely the determination of the f(a) func- 
tion and of the loss of information. Such consistency 
would be of great help for further understanding. These 
tasks remain, however, to be done. 

C. Mixing in many-body Hamiltonian systems 

Let us know address what might be considered as the 
Sancta Sanctorum of statistical mechanics, namely the 
systems that Boltzmann himself was mainly focusing on, 
i.e., classical many-body Hamiltonian systems. We know 
that no severe thermostatistical anomalies arrive for TV- 
body (i-dimensional systems whose (say two-body ) inter- 
actions are neither singular at the origin nor long-ranged. 
A most important example which violates both restric- 
tions is of course Newtonian gravitation. We shall here 
restrict to a simpler case, namely interactions that are 
well behaved at the origin, but which can be long-ranged. 
For specificity, let us assume a two-body potential which 
(attractively) decays at long distances like l/r Q with 
< a (for historical reasons we use a for this exponent; 
however, it has clearly nothing to do with the a appear- 
ing in the characterization of multifractals previously dis- 
cussed). The limit a — > oo corresponds to very short in- 
teractions; the other extreme (a = 0) corresponds to a 
typical Mean Field Approximation. No thermostatistical 
basic difficulties are expected for a > a c but the situation 
is quite different for < a < a c . As already mentioned, 
Fisher and collaborators established that a c = d for 
classical systems (a c > d is expected for quantum sys- 
tems). Therefore, the nontrivial case is < a/d < 1. 
Two physically different situations are to be understood, 
namely what happens at the limjv^oo lim^oo situation, 
and what happens at the lim^oo limjv->oo one. The two 
orderings are expected to be essentially equivalent for 
a/d > 1, and dramatically different otherwise. Let us 
first focus on the limjv—xjo lim^oo ordering. This should 
correspond to the traditional, BG-like, concept of ther- 
mal equilibrium, though with anomalous scalings. Let us 
be more specific and consider a fluid thermodynamically 
described by 

G(T, p, N) = U(T, p, N) - TS(T, p, N) + pV(T, p, N) , 

(11) 

hence, using the quantity N introduced in Section I, 



G{T,p,N) _ U{T,p,N) T S(T,p,N) | p V(T,p,N) 
NN NN N N N N 

(12) 

In the thermodynamic limit N — > oo we expect to obtain 
finite quantities as follows: 

g(f,p)=u(f,p)-fs(f,p)+pv(f,p), (13) 

with f = T/N and p = p/N. So, when a/d > 1, N is a 
constant, and we recover the usual concepts that G, U, S 
and V are extensive variables since they scale like N, 
whereas T and p are intensive variables since they are 
invariant with respect to N. If, however, < a/d < 1, 
the extensivity-intensivity concepts loose their usual sim- 
ple meanings. Indeed, there are now three, and not 
two, thermodynamical categories, namely the energetic 
ones (G, F, U, ...) which scale like NN, the pseudo- 
intensive ones (T, p, H, ...) which scale like N, and 
their conjugate pseudo-extensive variables (S, V, M, ...) 
which scale like N (if and only if expressed in terms of 
the pseudo- intensive variables). I can unfortunately not 
prove the general validity of the above scheme, but it has 
been shown to be true in all the models that have been 
numerically or analytically studied J5l[ (one- and two- 
dimensional Ising and Potts magnets, one-dimensional 
bond percolation and XY coupled rotators, one-, two- 
and three-dimensional Lennard- Jones-like fluids, among 
others). Also, for the case a = 0, we recover the tradi- 
tional scalings of Mean Field Approximations, where the 
Hamiltonian is artificially made extensive by dividing the 
coupling constant by N. By so doing, we are obliged to 
pay a conceptually high price, namely to have micro- 
scopic couplings which depend on N\ This long standing 
tradition in magnetism has, for very good reasons, never 
been adopted in astronomy: we are not aware of any as- 
tronomer dividing the universal gravitational constant G 
by any power of N in order to artificially make extensive a 
Hamiltonian which clearly is not! The thermodynamical 
scheme that has been described above escapes from such 
unpleasant criticism. I sometimes refer to the above dis- 
cussion of long-ranged systems as weak violation of BG 
statistical mechanics. Indeed, it becomes necessary to 
conveniently rescale the thermodynamic variables (using 
N) but the q — 1 approach remains the adequate one, 
i.e., energy distributions still are of the exponential type. 

Let us now address the much more complex 
lim^oo limjv^oo situation. At the present moment, the 
reliable informations that we have are mainly of numer- 
ical nature, more precisely computational molecular dy- 
namics. A paradigmatic system has been intensively ad- 
dressed during last years, namely a chain of coupled iner- 
tia! planar rotators (ferromagnetic XY-\ike interaction) , 
the interaction decaying as l/r a . The system being clas- 
sical, it is thermodynamically extensive for a > 1 and 
nonextensive for < a < 1. Important informations 
were established by Ruffo, in Florence, and collaborators 
p2[, who focused on the a = case. If we note En the 
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total energy associated with N rotators, it can be shown 
that a second order phase transition occurs at the crit- 
ical point e = E N /(NN) = 3/4 = e c . For e > e c the 
rotators are nearly decoupled and for e < e c they tend to 
clustcrize. Not surprising, e c depends on a, and vanishes 
above a critical value of a which is larger than unity. For 
a > 1, nothing very anomalous occurs. But, in the inter- 
val < a < 1, curious phenomena happen on both sides 
of e c (a). 

For e above e c (say e = 5), the (conveniently scaled) 
maximum Lyapunov exponent X max (N) goes to zero as 
N~ K , where k(o) was numerically established, by Celia 
Anteneodo and myself J53| (see Fig. 3). Nothing like 
this is observed for a > 1, where the numerical evi- 
dence strongly points that limjv-toc ^max(N) > 0. This 
of course suggests that, for a thermodynamically large 
system, mixing is not exponential, but possibly a power- 
law, i.e., weak mixing, as discussed above for the one- 
dimensional maps. The k(«) dependence has been re- 
cently studied by Andrea Giansanti, in Rome, and collab- 
orators for d = 2 and 3 as well. If ploted against a/d (see 
Fig. 3), an universal curve seems to emerge. This evidence 
is reinforced by the fact that it is possible to analytically 
prove, for this and other models, that k(0) = 1/3. It is 
consistently obtained k = for a/d> 1. 

Let us now focus on the anomalies below e c . The case 
e = 0.69 for a = has been and is being further studied 
in detail. The maximal Lyapunov exponent is here posi- 
tive in the limit N — * oo (for all < e < e c in fact) but, 
suggestively enough, the usual (BG) canonical ensemble 
equilibrium is preceded by a long, metastable-like, quasi 
stationary state (QSS), whose duration r diverges with 
AH This means that, if N is of the order of the Advo- 
gadro number, during a time that might be longer than 
the age of the universe, only this curious QSS will be ob- 
servable. So, we better understand it! As can be seen in 
Fig. 4, the average kinetic energy per particle e% depends 
on both N and time t. Only in the limjv^oo lim^oo the 
BG canonical temperature T is attained. If we take the 
limits the other way around a completely different result 
is obtained, namely approaches a finite value below 
T, and r diverges like N p . It becomes possible to spec- 
ulate that p{a) decreases from a finite positive value to 
zero, when a/d increases from zero to unity; thereafter, 
we of course expect p = 0, the system being then a well 
behaved BG one. To make this scenario more robust, 
consistent evidence is available concerning diffusion of 
the a — rotators. The numerical study of this phe- 
nomenon has shown J52|,^6| an average second moment 
of the angle increasing like t@ , where (3 depends on both 
N and t. For fixed N, (3 > 1 (i.e., superdiffusion) for a 
long time until, at t — {D stands for diffusion), makes 
a crossover to unity (normal diffusion). Analogously to 
what was described above, to increases like N a . Once 
again, it is possible to speculate that a (a) decreases from 
some finite positive value to zero when a/d increases from 
zero to unity. The door remains open to whether there 
is a simple connection between a and p. All these re- 



sults enable to conjecture that a scenario close to what 
is shown in Fig. 5 could indeed occur for some physical 
systems. Before closing this section, let me emphasize 
that what Fig. 4 strongly suggests is the possible validity 
of the zeroth principle of thermodynamics even out of a 
canonical BG scenario! Indeed, for large enough N, sys- 
tems can share the same " temperature" , thus being at 
thermal equilibrium, though this equilibrium is not the 
familiar one. This is quite remarkable if we take into 
account that there is nowadays quite strong numerical 
evidence that the distribution of velocities in such sys- 
tems is far from Gaussian, i.e., the familiar Maxwellian 
distribution of velocities is by all means violated. 



V. COMPARISON WITH EXPERIMENTAL 
RESULTS 

During this last decade many types of comparisons 
have been done between the present theory and experi- 
mental data. They have different epistemological status, 
and range from simple fittings (with sometimes clear, 
sometimes rather unclear physical insight), up to com- 
pletely closed theories, with no fitting parameters at all. 
Let me review a few among them. 



A. Turbulence 

X.-P. Huang and C.F. Driscoll exhibited in 1994 |§ 
some quite interesting nonneutral electronic plasma ex- 
periments done in a metallic cylinder in the presence of 
an axial magnetic field. They observed a turbulent ax- 
isymmetric metaequilibrium state, the electronic density 
radial distribution of which was measured. Its average 
(over typically 100 shots) monotonically decreased with 
the radial distance, disappearing at some radius sensibly 
smaller than the radius of the container (i.e., a cut-off 
was observed). They also proposed four phenomenolog- 
ical theories trying to reproduce the experimentally ob- 
served profile. The Restricted Minimum Enstrophy one 
provided a quite satisfactory first-approximation fitting. 
Bruce M. Boghosian, in Boston, showed in 1996 |5^] that 
the Huang and Driscoll successful attempt precisely cor- 
responds to the optimization of S q with q = 1/2, the 
necessary cut-off naturally emerging from the formalism. 
Since then, improved versions of his calculation p4] , p9| , 
as well as controversial arguments |jo), have been pub- 
lished. The merit of first connecting the present formal- 
ism to turbulence remains however doubtless. 

In recent months, Fernando M. Ramos and collabo- 
rators, in Sao Jose dos Campos, Brazil pH , Toshihiko 
Arimitsu and N. Arimitsu, in Tokyo |62|, as well as Chris- 
tian Beck, in London p3] |, have exhibited very interest- 
ing connections with fully developed turbulence. Beck's 
theory is a closed one, with no fitting parameters. He 
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obtained a quite impressive agreement with the experi- 
mental distribution of the differences (at distance r) of 
radial velocities, including its well known slight skew- 
ness (see Fig. 6). The entropic parameter is given by 
l/(q — l) = l+log 2 (r/?7) where 77 is the viscosity. The fact 
that q is not universal certainly reminds us the critical- 
ity of the two-dimensional XY ferromagnet, with which 
nonextensivity presents in fact various analogies. 

B. High-energy collisions 

Ignacio Bediaga and collaborators (64|], in Rio de 
Janeiro, have worked out, along Rolf Hagedorn's 1965 
lines |65| , a phenomenological theory for the distribution 
of transverse momenta in the hadronic jets emerging from 
electron-positron annihilation after central collisions at 
energies ranging from 14 to 161 Gev. Indeed, since early 
ideas of Fermi, and later of Field and Feynman, a ther- 
modynamical equilibrium scenario has been advanced for 
this distribution of transverse momenta. Hagedorn de- 
veloped a full theory based on the BG distribution. The 
success was only partial. Indeed, a central ingredient 
of the physical picture is that higher collision energies 
do not increase the transverse momenta temperature T 
but increase instead the number of involved bosons that 
are produced (like water boiling at higher flux of energy, 
where only the rate of vapor production is increased, 
but not the temperature). It is this central physical ex- 
pectance that was not fulfilled by Hagedorn's calculation; 
indeed, to fit the curves associated with increasingly high 
energies, increasingly high values for T become necessary. 
Furthermore, an inflexion point clearly emerges in the 
distribution, which by no means can be reproduced by 
the BG exponential distribution. The adequation of only 
two parameters (T and q) have enabled Bediaga et al |64|] 
to fit amazingly large sets of experimental data (see Fig. 
7). And, in order to make the expected picture com- 
plete, the central demand of T being independent from 
the collision energy is fulfilled! 

Grzegorz Wilk, in Warsaw, and collaborator s [p6[ , as 
well as D.B. Walton and J. Rafelski in Tucson J67|, have 
provided further evidences of the applicability of the 
present ideas to high energy physics, where nonmarko- 
vian processes and long-range interactions are commonly 
accepted hypothesis. Consistent evidence has also been 
provided by D.B. Ion and M.L.D. Ion, in Bucharest, very 
especially in hadronic scattering H] . 

C. Solar neutrino problem 

Piero Quarati and collaborators, in Torino and 
Cagliari, have been advancing since 1996 |39| an inter- 
esting possibility concerning the famous solar neutrino 
problem. Indeed, it is known that calculations within the 
so-called Standard Solar Model (SSM) provide a neutrino 



flux to be detected on the surface of the Earth which is 
roughly the double of what is actually measured in several 
laboratories around the world. This enigmatic discrep- 
ancy intrigues the specialists since over 20-30 years ago. 
At least two, non exclusive, possible causes are under cur- 
rent analysis. One of them concerns the nature of neu- 
trinos: they could oscillate in such a way that only part 
of them would be detectable by the present experimental 
devices. The second one addresses the possibility that the 
SSM could need to be sensibly improved in what concerns 
the production of neutrinos at the solar core plasma. It is 
along this line that Quarati's suggestion progresses. In- 
deed, the neutrino flux is related to the total area of the 
so-called Gamow peak, which is in turn due to the prod- 
uct of the thermal equilibrium BG distribution function 
(which decreases with energy) and the penetration factor 
(which increases with energy) . The position of the peak 
is at energies 10 times larger than kT, therefore only the 
far tail of the thermal distribution is concerned. Quarati 
et al argue that very slight departures from q = 1 (of the 
order of 0.01) are enough to substantially modify (close 
to the desired factor two) the area of the Gamow peak. 
They also show that this degree of departure remains 
within the experimental accuracy of other independent 
measurements, such as the helioseismographic ones. This 
slight nonextensivity would be related to quite plausible 
nonmarkovian processes and other anomalies. 

D. Others 

Many other phenomena have been, during recent years, 
temptatively connected to the present nonextensivity. 
Let me briefly mention some of them. The experimental, 
clearly non Gaussian, distribution of velocities in various 
systems has been shown to be very satisfactorily fitted 
by q ^ 1 curves. Such is the case of the distribution 
(observed from the COBE satellite) of peculiar velocities 
of spiral galaxies, shown f70fl , in 1998 by Quarati, myself 
and collaborators, to be very well fitted by q ~ 0.24. Such 
is also the case of the distribution that James Glazier, in 
Notre Dame University, and collaborators [JFlJ , measured 
for the horizontal velocities of Hydra Vulgaris living in 
physiological solution. They were remarkably well fitted 
using q ~ 1.5. Such is finally the case of the computa- 
tionally simulated granular matter involving mesoscopic 
inelastic collisions. Indeed, the distribution function that 
Y.-H. Taguchi and H. Takayasu used to fit their 1995 data 
p"2| precisely is a typical q > 1 one. 

In 1995 |7^| , F.C. Sa Barreto, in Belo Horizonte, Brazil, 
E.D. Loh, in Michigan State University, and myself ad- 
vanced the possibility of explaining very tiny departures 
from the black-body radiation Planck's law for fitting 
the COBE satellite data for the cosmic background mi- 
crowave radiation. Indeed, values of q departing from 
unity by the order of 10 -5 were advanced. This pos- 
sibility has since then be further analyzed by a variety 
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of scientists, the conclusion still remaining basically the 
same. Only more precise (perhaps 10 times more precise 
would be enough) experimental data would confirm or 
reject such a possibility, physically motivated by present 
or ancient subtle effects of gravitation, which could even 
lead ultimately to a modification of our understanding 
of the nature of spacetime (at very small scales, it could 
very well be noncontinuous and even nondifferentiable, 
in contrast with our usual perceptions!). 

Molecules like CO and O2 in some hemoproteins can be 
dissociated from their natural positions by light flashes, 
as many classical experiments have shown. Under specific 
circumstances, these molecules then tend to re-associate 
with a time-dependent rate. The experimental rate has 
been shown, in 1999 by G. Bemski, R.S. Mendes and 
myself |73|] , to be very well fitted by the solutions of 
a nonextensive-inspired differential equation. The in- 
trinsic fractality of such proteins is believed to be the 
physical motivation for such approach. Numbers of cita- 
tions of scientific publications can also be well fitted with 
q > 1 curves, consistently with S. Denisov's 1997 recov- 
ering, within the present formalism, of the so-called Zipf- 
Mandelbrot law for linguistics j74|. Oscar Sotolongo- 
Costa, in La Habana, and collaborators jf5|, as well as 
K.K. Gudima, in Caen, France, and collaborators [ [76| , 
have recently invoked the present formalism for describ- 
ing multifragmentation. 

Last but not least, it is worthy mentioning that the 
present framework has induced [ J77[ quite quicker versions 
of the Simulated Annealing methods for global optimiza- 
tion. This methodology has found successful applications 
in theoretical chemistry, in particular for studying pro- 
teins, their folding and related phenomena. Such pro- 
cedures are currently being used by John E. Straub, at 
the Boston University, and collaborators, by Kleber C. 
Mundim, in Salvador, Brazil, Donald E. Ellis, in Chicago, 
and collaborators, by Yuko Okamoto, in Okazaki, Ulrich 
H.E. Hansmann, in Houghton, and collaborators, among 
others @. 

VI. FINAL COMMENTS 

Many issues remain, at the present moment, partially 
or fully open to better understanding concerning the 
present attempt of adequately generalizing Boltzmann- 
Gibbs statistical mechanics and thermodynamics. These 
include (i) the direct checking of the energy distribu- 
tion which generalizes the BG factor (and consistently 
the connection between q and (a, d) for long-range inter- 
acting Hamiltonian systems (both finite and large), and 
similar connections for nonconservative ones); (ii) appro- 
priate understanding within the renormalization group 
framework (Alberto Robledo, in Mexico, has recently ini- 
tiated this line frSj ]; see also pO[); (iii) possible firm con- 
nections of this formalism with the quantum group for- 
malism for generalizing quantum mechanics (the great 



analogies that these two nonextensive formalisms exhibit 
are since long being explored, but not yet deeply under- 
stood); (iv) possible connections with the Levy- like prob- 
lematic reviewed recently M] by George M. Zaslavsky; 
(v) the strict conditions (in terms of the N — > 00 and 
t — > 00 limits) under which the zero'^-principle of ther- 
modynamics (i.e., the criterium for thermal equilibrium) 
holds; (vi) the microscopic and mathematical interpreta- 
tion of the direct and escort distributions (which possibly 
reflect the choice of a (multi)ffactal description, or of a 
description in terms of Lebesgue-integrable quantities); 
(vii) the general connection between symmetry and en- 
tropic form adapted to measure information about sys- 
tems having that particular symmetry; (viii) the rigorous 
basis for the generalization herein described of concepts 
such as Kolmogorov-Sinai entropy, Lyapunov exponents, 
Pesin equality, escape rates, etc; (ix) the possibility of 
having, for long-range interacting systems, a special value 
for q which, through the use of appropriately rescaled 
variables, would enable a new kind of extensivity in the 
sense that S q (A + B) = S q {A) + S q (B), where the large 
subsystems A and B would of course be not independent 
(in the sense of theory of probabilities) but, on the con- 
trary, strongly coupled; (x) the connection between the 
index q appearing in properties such as the sensitivity 
to the initial conditions, and that appearing in proper- 
ties such as the microscopic distribution of energies (the 
most typical values for the former appear to be smaller 
than unity, whereas they are larger than unity for the lat- 
ter) ; (xi) the direct checking of the type of mixing which 
occurs in the phase-space of many-body long-range in- 
teracting large systems; (xii) the deep connection with 
quantum entanglement through its intrinsic nonlocality. 
The last point certainly appears as very promising and 
relevant. An important step forward was very recently 
done by Abe and Rajagopal Q; indeed, they succeeded 
in recovering, through use of the present formalism, the 
x = 1/3 Peres (necessary and sufficient) condition, which 
is known to be the strongest one (stronger than Bell in- 
equality) for having local realism (some degree of sepa- 
rability of the full density matrix into those associated 
with subsystems A and B) for the bipartite spin-i sys- 
tem (Werner-Popescu state). The probability of this ap- 
proach being fundamentally correct is enhanced nowa- 
days by contributions produced from quite different view- 
points Q . 

As we can see, the number of important issues still 
needing enlightening is quite large. This makes the en- 
terprise but more stimulating! However, caution is rec- 
ommended in spite of the sensible amount of objective 
successes of the proposal. Indeed, given the enormous 
impact of the thermodynamical concepts in Physics, one 
should address all these points with circumspection. Nev- 
ertheless, an important new perspective seems to be al- 
ready acquired. Nikolai S. Krylov showed, half a cen- 
tury ago, that the key concept in the foundation of stan- 
dard statistical mechanics is not ergodicity but mixing. 
Indeed, he showed that the Lyapunov exponents essen- 
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tially control the relaxation times towards BG equilib- 
rium. The entropy associated with such systems is vastly 
known to be the BG logarithmic one. What emerges now 
is an even more fascinating possibility, namely the type 
of mixing seems to determine the microscopic form of the 
entropy to be used. If the mixing is of the exponential 
type (strong chaos), then the entropy of course is the BG 
one (i.e., q = 1). If the mixing is of the power type (weak 
chaos), then the mixing exponent (1/(1 — g) apparently) 
would determine the anomalous value of q to be used for 
the entropy, which would in turn determine the thermal 
equilibrium distribution, as well as all of its thermody- 
namical consequences! The kingdom of the exponentials 
would then be occasionally replaced by the kingdom of 
the power laws, with their relevance in biology, economics 
and other complex systems! 
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APPENDIX 

Scene at the restaurant 

A U LABYRINTHE DES ENTROPIES 

- Would you have some fresh entropies today, for me 
and my friends ? 

- Absolutely Sir ! We have them extensive or not, with 
definite concavity or not, nonnegative defined or other- 
wise, quantum, classical, relative, cross or mutual, in- 
cluded in several others with a small supplement, com- 
posable or not, expansible or not, totally optimized or a 
little rare, even completely out of equilibrium... single- 
trajectory-based or ensemble-based... 

You can have them at the good old Boltzmann mag- 
nificent style, doree a la Gibbs, very subtle, or von Neu- 
mann..., with pepper a la Jaynes, or the popular Shan- 
non, oh, my God, I was forgetting the esoteric, superb, 
macroscopic Clausius, the surprising Fisher, the refined 
Kolmogorov-Sinai, a la Kullback-Leibler for comparison, 
Renyi with multifractal dressing, with cybernetic sauce a 
la Harvda and Charvat, Vajda and Daroczy, or even the 
all-taste Sharma and Mittal... 

In 1988 we started serving them with Brazilian touch, 
if you wish to try, it leaves a tropical arriere-gout in your 
mouth! And since then, our chefs have introduced not less 
than ten new recipes... Curado, with exponentials, An- 
teneodo with tango flavor, Plastino 's, excellent as fam- 
ily dish, Landsberg, Papa, Johal with curry, Borges and 



Roditi, Rajagopal and Abe... 

They are all delicious ! How many do I serve you to- 
day? 

FIGURES 

Fig. 1 - Equilibrium-like probability versus energy 
for typical values of q, namely, from top to bottom at 
low energies, q = 0, 1/4, 1/2, 2/3, 1, 3, oo (the latter 
collapses onto the ordinate at the origin, and vanishes 
for all positive energies). For q = 1 we have Boltzmann 
exponential factor; for q > 1 we have a power-law tail; 
for q < 1 there is a cutoff above which the probability 
vanishes. 

Fig. 2 - Time evolution of S q : (a) a = 2 (q* = 1); (b) 
a = a c (q* = 0.2445; R measures the degree of nonlin- 
earity of the curves in the relevant intermediate region). 

Fig. 3 - Coupled planar rotators: k versus a/d (d= 1: 
f33| |; d — 2,3: pjj). The solid line is a guide to the eye. 

Fig. 4 - Coupled planar rotators (a = 0, d = 1): Time 
evolution of T. From Latora and Rapisarda p5|). 

Fig. 5 - Conjectural time dependence of the prob- 
ability p of having an energy E for classical systems 
(N = N* + 1). 

Fig. 6 - Distribution of radial velocity differences in 
fully developed turbulence (from top to bottom: r/rj 
3.3, 23.6, 100). From Beck " 



Fig. 7 - Distribution of hadronic transverse momenta: 
(a) cross-section (dotted line: Hagedorn's theory; solid 
lines: g ^ 1 theory )j_(b) Fitting parameters q and To. 
From Bediaga et al |34| . 
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